Foundations of Physics Letters, Vol. 6, No. 6, p. 571-590 (1993) 



APPARENT WAVE FUNCTION COLLAPSE 
CAUSED BY SCATTERING 



Max Tegmark 

Department of Physics, University of California 
Berkeley, California 94720 



Received April 19, 1993; revised June 10, 1993 



Some experimental implications of the recent progress on wave function col- 
lapse are calculated. Exact results are derived for the center-of-mass wave 
function collapse caused by random scatterings and applied to a range of 
specific examples. The results show that recently proposed experiments 
to measure the GRW effect are likely to fail, since the effect of naturally 
occurring scatterings is of the same form as the GRW effect but gener- 
ally much stronger. The same goes for attempts to measure the collapse 
caused by quantum gravity as suggested by Hawking and others. The re- 
sults also indicate that macroscopic systems tend to be found in states with 
AxAp = but microscopic systems in highly tiltedly squeezed states 

with AxAp h. 
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1. INTRODUCTION 

The problem of how to interpret measurement in quantum mechanics 
has caused intense debate ever since 1926 and shows little sign of abating. 
A whole slew of interpretations have been proposed and can be divided into 
two main categories: collapse theories and non-collapse theories. 

In the former category, one of the most successful to date is that pro- 
posed by Ghirardi, Rimini and Weber (GRW) in 1986 [1], which shows that 
both micro- and macroscopic systems can be described by the same dy- 
namical equation providing that an extra term is added to the Heisenberg 
equation of motion for the density matrix: 



p = -jr[H,p]-A(p-T\p]) 



(1) 



They show that if this ad hoc third term is added and if T[p] is chosen 
in a particular way that singles out the position representation as special, 
then the usual problems regarding superpositions of macroscopic systems 
disappear if A and a second parameter are chosen appropriately. This theory 
has subsequently been generalized. 

In the second category of theories, perhaps the most radical is the one 
proposed by Everett, Wheeler, Cooper, DeWitt, and others [2-7] between 
1957 and 1970, which shows that even if one assumes that the wave func- 
tion containing the observer evolves causally according to the Schrodinger 
equation, the observer will subjectively experience wave function collapse. 
Zeh, Kiibler, Joos, Machida, Namiki, Zurek, Unruh, Cini, Peres, Partovi, 
Gallis, Fleming, Hartle and others [8-20] have strengthened this position 
by showing that for macroscopic objects, their inevitable interaction with 
the environment leads to a dynamic reduction of the density matrix (what 
is widely known as wave function collapse) and superselection rules. These 
superselection rules tend to favor "classical" states, and explain why we 
never experience say spatial superpositions of cars or superpositions of liv- 
ing and dead cats. Hence this interaction with the environment shows why 
the position operator and its eigenstates play such an important role in 
our perception of the world, even though the position operator is a priori 
merely one out of a family of infinitely many self-adjoint operators. This 
will be referred to as the decoherence effect. 

Unfortunately, rather scant attention has yet been given to the ex- 
perimental implications of the decoherence effect and to actual physical 
parameters. This paper addresses such practical issues, focusing on scat- 
tering. Recently, experiments have been proposed [21, 22] to try to detect 
the GRW effect, but this paper shows that such experiments are likely to 
fail, since a GRW effect with the parameters originally proposed would be 
entirely drowned out by environmental noise. 

In Sec. 2, the effect of a single scattering is calculated and shown to 
damp the off-diagonal elements in the reduced spatial density matrix by a 
factor that is simply the Fourier transform of the probability distribution 
for different momentum transfers. In Sec. 3, the Heisenberg and Wigner 
equations of motion are modified to incorporate these usually neglected ef- 
fects. In Sec. 4 the results are applied to a variety of cases of physical 
interest and compared to the predictions of the GRW theory and quantum 
gravity. It is seen that scattering and the GRW effect have almost identi- 
cal effects on the reduced density matrix, although the interpretations are 
completely different. Finally, Sec. 5 contains a brief discussion of what 
interpretational problems of quantum mechanics the decoherence approach 
does and does not solve. 



2. THE EFFECT OF A SINGLE SCATTERING 

As an introduction to the calculations in this section, consider the fol- 
lowing simple example of decoherence: A spin- ^ silver atom is prepared with 



its spin in the x-direction, and then somebody measures its spin component 
in the z direction without telling us the result. This changes our density 
matrix for the atomic spin from pi to pf, where in the z-representation 



Thus the density matrix is reduced from describing a pure state to a mixed 
state. This prediction is common to all interpretations of quantum mechan- 
ics, but arrived at in two conceptually very different ways: Collapse-theories 
postulate that the time-evolution of the wave function of the universe is not 
governed by the Schrodinger equation during measurement, but changes dis- 
continuously and non-causally so that afterwards the spin really is up or 
down in the z-direction - we just do not know which. Non-collapse theories 
compute the same density matrix pf by letting the total system of observer 
and observed evolve according to the Schrodinger equation with a Hamil- 
tonian such that they become perfectly correlated, and take a partial trace 
over the observer degrees of freedom to obtain pf. These two incompatible 
viewpoints are often referred to as Heisenberg reduction and Von Neumann 
reduction, respectively. (A detailed discussion of these matters is given by 
Everett [3] and Kraus [23].) Here we will adopt the latter approach, and 
refer to it as the decoherence approach. 

Although it is often convenient to treat particles as isolated systems, 
we all know that this is merely an approximation. Occasionally a photon 
from the sun scatters off of our "isolated system" . More difficult to shield 
experiments from are muons created by cosmic rays, cosmic neutrinos, the 
300K blackbody radiation from our surrounding and radiation from traces 
of radioactive isotopes in the materials that our measurement apparatus is 
made of. All these events change the density matrix of our particle. Joos 
and Zeh have studied such scattering effects in the macroscopic limit by 
ignoring recoil [11], whereas the following treatment applies also to micro- 
scopic systems. 

Let us chose as our system a nonrelativistic particle of mass m whose 
location is described by a density matrix p and whose time-evolution would 
be governed by a Hamiltonian H if it were truly isolated from its environ- 
ment. Let the inevitable interaction with the environment be given by an 
interaction Hamiltonian Hint- In this paper we will limit our attention to 
the special class of interactions with the environment that can be treated as 
isolated scattering processes. By this we mean that Hint 7^ only during 
time intervals much shorter than the dynamical timescale of the system we 
are studying, so that we can approximate the change in p as instantaneous 
and given by a transition matrix T, 



where p T is the density matrix for the total system of our particle and 
an external particle that scatters off of it. This is normally a good ap- 
proximation when our system interacts with a rapidly moving particle in 




and 




pf - p T f := TpjT\ 



its vicinity. For instance, it takes a photon only about 10 seconds to 
traverse an atom. 

We will make the following assumptions about the T-matrix and the 
initial data: 

Assumption (I) T conserves energy and momentum. (This is equiv- 
alent to T being invariant under temporal and spatial translations.) Let 
|pk) denote the state where our system has momentum p and the incident 
particle has momentum k. Then (I) implies that 

(p'k'l T |pk> = 5(p' + k' - p - k)a pk ( P ' - p), (2) 

where a p k(q) is the probability amplitude for the momentum transfer to our 
system to be q. This function is independent of time by energy conservation 
but may depend on both p and k. 

Assumption (II) The function a p k is independent of p, i.e of the 
motion of our system. Hence we will write it as a^. (This is a good ap- 
proximation if the velocity of the incident particle is much greater than the 
velocity spread in p). 

Assumption (III) The incident particle is in a momentum eigenstate 
or an incoherent mixture of momentum eigenstates. (The linewidth must 
be much smaller than the wavelength. For the photons we observe, the 
linewidth is typically less than 1% of the wavelength.) 

To avoid normalization problems, let us first restrict ourselves to L 2 
functions with periodic boundary conditions on a cube of volume V. Unless 
otherwise specified, all integrals below are to be taken over this cube and 
all sums are to be taken over the discrete set of vectors 

ft := {(2ir/V 1/3 )(n x ,n y ,n z )\n x ,ny,n z integers}. 

To conserve probability, Ok must be normalized so that -Pk(q) : = |°k(q)| 2 
is a probability distribution over q, i.e. 

q 

(In the literature, T is often normalized so that this sum equals the total 
cross section a instead. We will take account of the cross section in Sec. 3.) 

Equation (2) shows that in the position representation, T transforms 
the state |pk) from e l P- x i e lk - x 2 m to 

5^«k(q)e i(p+q) - Xl e i(k - q) - X2 = e ip - Xl e ik - X2 a k (x 2 - Xl ), 
q 



where is the discrete Fourier transform of ak, so for a normalized ini- 
tial two-particle wave function V>i(xi,x 2 ) = </>(xi)y _1 / 2 e lk ' X2 , we have by 
linearity that 

r^(x 1 ,x 2 ) = r(</»(x 1 )y- 1 / 2 e ik - X2 ) = 0(x 1 )y- 1 / 2 e ik - X2 a k (x 2 - Xl ) 
= ^i(xi,x 2 )a k (x 2 - xi) 

for any one-particle wave function </>(xi). Let us use the notation p(x, y) := 
(y| p |x) for density matrices and units where ft = 1. 

Theorem p/(x, y) = pi(x, y)Pk(y — x), where Pk is the Fourier trans- 
form of Pk := |«k| • 

Proof The reduced density matrix p for our particle is obtained by 
taking a partial trace of the density matrix p T of the two-particle system, 
so 

Pi(x,y) = |v*(x,z)^(y,z)d 3 z = |^(x)0(y)id 3 z = ^(x)0(y), 
and by the above, 

p/(x,y) = y ^/( x > z )^/(y,z)d 3 ^ 

(^i(x,z)a k (z - x))*(^i(y,z)a k (z - y))d 3 z 
= ^ (x)0(y)-^a k (z-x)a k (z-y)d 3 z 

= Pi(x,y)^ y a k ( z -x)a k (z-y)(i 3 2;. 

The not very elegant V made its last appearance in this paper on the previ- 
ous line. Letting it approach infinity and working with Fourier transforms 
instead of Fourier series from here on, the last equation turns into 

P/(x,y) =Pi(x,y)^Tg y a*(z-x)a k (z-y)(f 3 z. 

By substituting u = z — y and using d k (u) = a k (— u), we obtain 
P/(x,y) = pi(x,y)^3 y a*(y-x-u)a k (u)d 3 u 

= P;(x,y)^^(a k *a k )(y-x). 



/ 



Using the convolution theorem / * g = (2iT) 3 fg now yields 



p/(x,y) = p l (x,y)^-^(2vr) 3 a*a k (y-x) = p;(x, y)P k (y - x), 

which completes the proof for the case where our system is initially in a pure 
state, i.e., where p^(x, y) can be written in the form <f>* (x)0(y). Since an 
arbitrary density matrix can be written as a sum of density matrices for pure 
states, the proof for the general case follows directly from superposition. 

Corollary If the incident photon is described by a density matrix 
diagonal in the momentum representation, say an incoherent superposition 
of plane waves with the momentum probability distribution given by ^(k), 
then 

p,(x,y) = p i (x,y)P(y-x), where P(q) := J P k (<i)p(k)d 3 k. (3) 

This again follows directly from superposition. 

Thus we see that the net result of this interaction of our "isolated sys- 
tem" with the outside world is simply to multiply its spatial density matrix 
by a function. Note that this function does not depend on the complex 
amplitude Ok itself, but only on its squared modulus, the probability distri- 
bution for different momentum transfers. The latter is uniquely determined 
by the differential scattering cross section a(9, <p), so a is the only physical 
input we will need to calculate how p evolves over time. 

Before we turn to calculating P for specific physical examples, let us 
make a few observations about Eq. (3) that are valid for an arbitrary 
probability distribution P. It is straightforward to prove the following: 

Observation (I) |P(x)| < 1, with equality if x = or 
P(q) = <5(q-q ). 

Observation (II) tr pf = tr p. 

Observation (III) tr p 2 , < trpf, with equality iff pi is diagonal or 

P(q)=<y(q-qo). 

Observation (IV) P(x) — > as |x| — > oo if P is an integrable function. 

Observation (V) If we define Q to be the mean and S to be the 
covariance matrix of the probability distribution P, then 

P(x) 1 iQm%m 2^Q m Q n ^mn)X m X n -\- 0(|x| ). 

(Repeated indices are to be summed over, from 1 to 3.) 

Observation (I) tells us that all off-diagonal elements of the density 
matrix will be damped if there is any uncertainty in the outcome of the 
scattering. Observation (II) simply shows that probability is conserved, 
whereas (III) says that the density matrix generally becomes less pure - we 
recall that trp 2 = 1 for pure states, whereas trp 2 takes its (nonnegative) 



minimum value for states of which we have zero knowledge. (IV), which 
is known as Riemann-Lebesgue's Lemma, tells us that for most physically 
realistic cases, the density matrix elements very far from the diagonal get 
almost entirely damped out. (V) gives us a good grip on how the density 
matrix changes near the diagonal, which will prove useful in the following 
section. 



3. THE MODIFIED HEISENBERG AND WIGNER EQUA- 
TIONS OF MOTION 

Above we derived the effects of a single scattering event. Now we 
will show that exposure to a constant particle flux amounts to a simple 
modification of the Heisenberg equation of motion for our density matrix p. 

Let a denote the total scattering cross section and <E> the average par- 
ticle flux per unit area per unit time. We know from experiment that 
the temporal distribution of scattering events is well modeled by a Poisson 
process with intensity A := Since the probability of one scattering 

occurring during the infinitesimal time interval dt is Adt and that for no 
scattering is 1 — Adt, by Eq. (3) we would get 

p(x, y, t + dt) = p(x, y, t)P(y - x) Adt + p(x, y, t) (1 - Adt), 

i.e. p(x,y,t) = — A(l — P(y — x))p(x, y, t) if our scattering were the only 
process that changed p. Since p is also changed by its normal Sshrodinger 
time evolution, we obtain the following master equation: 

p = -±[H,p]-A(p-T[p]), (4) 

where 

(y|T[p]|x) :=P(y-x)p(x,y). 

We notice that this is of the same form as the GRW equation (1), and soon 
we will indeed see that the T[p] that we have derived is quite similar that 
postulated by GRW. They postulate that (y| T[p] |x) is a Gaussian with 

standard deviation Ao- We will refer to (1 — P) as the decoherence func- 
tion (what Gallis and Fleming call the decorrelation factor) . This function 
summarizes all there is to know about the scattering environment, since 
(4) shows that it together with the Hamiltonian specifies the dynamical 
behavior of our system entirely. 

Let us take a closer look at the last term. In the position representation, 
observation (IV) shows that it approaches — Ap(x, y) far from the diagonal. 
Near the diagonal, we can use (V) to expand it as 

(y|-A(p-T[p]) |x) « -iAQ-(y-x.)-±A(Q m Q n +S mn )(y m -x m )(y n -x n ), 



where we have dropped cubic and higher order terms in |y — x|, the distance 
to the diagonal. Now since (y| [Q • x, p] |x) = Q • (y — x)p, we can absorb 
the first term into the Hamiltonian as a linear potential and write (4) as 



where < y|-D[p]|x \{Q m Q n + S mn )(y m - x m )(y n - z n )p(x,y). 

This linear potential should come as no surprise — it is simply the radiation 
pressure term, and causes no dissipation. (By dissipation we will mean a 
process that increases entropy, i.e., that converts pure states into mixed 
states. More quantitatively, we will say that we have dissipation if the 
linear entropy 1 — tr p 2 increases.) From observation (III) we know that in 
general we do have dissipation, so the cause of this must be the term D[p\. 

To get a better feeling for what is happening, let us look at what effect 
a single scattering would have in phase space. By transforming Eq. (3) to 
the Wigner representation [24, 25] 



and doing some algebra, we see that a single scattering has the effect 



i.e. a smearing out in momentum space. Thinking of W as a probability 
distribution in phase space, this convolution with the probability distribu- 
tion P(q) corresponds to giving our particle a random momentum kick. 
Thus we can use the central limit theorem and approximate the effect of 
n consecutive hits by an equation identical to (5) but with P(q) replaced 
by a Gaussian with mean nQ and covariance matrix nS. Since the exact 
number of scatterings is not known but Poisson distributed with mean and 
variance equal tow = At, the Gaussian will in fact have the covariance ma- 
trix At(Q m Q n + S mn ). Thus we see that for the non-isotropic case Q^O 
the purely epistemological uncertainty as to how many scatterings have oc- 
curred increases the rate of wave function collapse, the rate of damping of 
off-diagonal elements. 

Returning to our density matrices, this means that if A is so large that 
there are many scatterings on time scales shorter than that of ordinary 
Schrodinger evolution, then we can replace the function P in the master 
equation (4) by a Gaussian with the same mean and covariance matrix. 
Equivalently, we can replace P by a Gaussian with the same values of ze- 
roth, first and second derivatives at the origin. This shows that in the limit 
of large A, our Eq. (4) reproduces the GRW equation (1) exactly. Further- 
more, since the Gaussian is the Green function of the diffusion equation, 



P = 



UH + AQ- x ,p]-AD[p], 






this means that we can incorporate our dissipative term into the Wigner 
equation [24, 25] as a simple diffusion term of the type V 2 : 



where V(x) is the potential and the diffusion coefficient := A(5 , ij - + 
QiQj). Here the spatial derivatives in the infinite sum are to be understood 
to act only on V, not on W. There are two well-known limits in which the 
Wigner equation goes over into the Liouville equation of classical statistical 
mechanics: when V is at most quadratic, and when h — > 0. Because of 
our extra diffusion term, we get yet a third classical limit: in the limit 
of large D^, the diffusive smoothing becomes so effective that it damps 
out all the momentum-derivatives in the infinite sum, and (6) approaches 
the Liouville equation with diffusion, an equation of Fokker-Planck type. 
This is yet another example of how macroscopic objects start behaving 
classically, since as we will soon see, Djj is roughly proportional to the size 
of our object. Thus an object will evolve according to classical dynamics 
if it has a strong interaction with its environment. (The diffusion term 
in the resulting Liouville equation is in no way a departure from classical 
dynamics, since the Brownian motion due to random scatterings must be 
taken into account also in a purely classically analysis.) 

4. SPECIFIC EXAMPLES 

In this section we will apply our results to specific scattering processes. 
We will first calculate the shape and width of the decoherence function, then 
study what the density matrix of a free particle converges to as t — > oo. 

4.1. The Decoherence Function 

Now let us evaluate Q, Dij and the decoherence function 1 — P 
for some physically interesting cases. If the scattering cross section is given 
by da/dQ, = a f(6,ip), where a is the total cross section and the angular 
part / is normalized so as to integrate to unity, then -Pk(q) = <5(|q — k| — 
k) f /k 2 , where k := |k|. Using the properties of the Fourier transform under 
translation and reflection and choosing our coordinates so that r = (0, 0, r), 
we get 








With a generic anisotropic radiation spectrum //(k), the rate of decoher- 
ence will be different along different spatial directions and we will get ra- 
diation pressure. Since neither of these two complications is particularly 
illuminating, we will restrict ourselves to isotropic radiation, i.e. take 
/x(k) = (47rk 2 ) _1 Ao^(Ao|k|), where the spectrum v is a probability distribu- 
tion on the positive real line and Ao is some typical wavelength. Performing 
the angular integration, this yields 

/■oo 

P(r) = / g(ur/X )u(u)du, (7) 
J o 

sinx f 1 

g(x) := / e txu f(&r ccos u,(p)dud(p. 

x Jo 

For the case of photon scattering against both a free charge and a dielectric 
sphere much smaller than the photon wavelength, we get [26] the angular 
dependence 

which yields 

q r 

cosx + (x 2 — 1)- 



where 



3 

g(x) =gi(x) := - 



x 



sin 2 x 



x 4 



(8) 



We will refer to this mixture of S-waves and D-waves as SD-wave scattering. 
Another physically important case is pure S-wave scattering, i.e. 

1 { sin x \ ^ 

Z( 6 '' V 3 ) = ^ which yields g(x) = g 2 (x) := I 1 . (9) 

This case applies among other things to an opaque spherical object of ra- 
dius a much larger than A e //, which we can use to model say a dust particle 
scattering optical photons. Here the total cross section is frequency inde- 
pendent, roughly equal to the geometrical cross section ira 2 , and perhaps 
surprisingly, the scattering amplitude turns out to be the same in all direc- 
tions. 

Since we are restricting ourselves to isotropic radiation, we simply have 
the mean Q = and the covariance matrix is proportional to the identity 
matrix, i.e. Sij = s 2 5ij, so all we need to calculate is the standard deviation 
s. From Eq. (7) we get 

/>oo 

P"(0) = A " 2 / v(x)x 2 g"(0)dx = g"(0) (x 2 ) , 



so 



s = X^[-g"(0)]^ 2 (x 2 ) 1 



(10) 



since the variance of a probability distribution of zero mean is the negative of 
the second derivative of its Fourier transform at the origin. For the functions 
in (8) and (9), the values we need are -<(0) = 11/15 and -g%(0) = 2/3. 
Thus apart from these numerical parameters depending only on the angular 
part of scattering cross section, we see that the standard deviation of the 
momentum kick is simply a certain spectrally averaged momentum. 
Let us define the effective wavelength as 

A e// := l/s = \ [-g"(OT^(x>)- 1/2 . 

Then the diffusion matrix Djj = A5ij , where the scalar diffusion coefficient 
A := A/Agjj. Sunlight on earth, the 300-ftT radiation from our surrounding, 
the cosmic microwave background radiation and the cosmic neutrino back- 
ground all have Planck spectra, corresponding to temperatures of roughly 
5800-fT, 300-ftT , 2.7K and 2. OK, respectively. For a Planck spectrum, we 
have 

2 X 2 

u ( x ) = cm e x _ i ' 

where £(s) := Ylri=i n ~ s IS the Piemann Zeta function. For this particular 
case, (x 2 ) = 4!£(5)/2!((3), so we get A e // ~ 0.381Ao for the S-wave case and 
X e ff w 0.363Ao for the SD-wave case. For a point spectrum u(x) = 5{x— 1) 
with a single wavelength Ao we simply get (x 2 ) = 1, so the corresponding 
values are A e // ps 1.225Ao and A e // « 1.168Ao- 

Now let us calculate the decoherence function 1 — P for the Planck 
case. By expanding v as a geometric series, we get 

P(r) = / g(ur/\ )u{u)du = — — V / g{xu)u 2 e- nu du. 
Jo ^) n=1 Jo 

For the S-wave case g = g2, this integral can be done elementarily, yielding 




Numerical integration is required in the case g = gi, and gives a decoher- 
ence function 1 — P(r/A e //) that differs by less than 1% from the S-wave 
case when appropriately rescaled. This decoherence function is plotted for 
Figure 1, together with the S-wave and SD-wave decoherence functions for 
point spectra. Also plotted is the GRW decoherence function, for which P 
is Gaussian. 

For dielectric spheres with a <C X e ff and frequencies well below any 
resonances, the cross section depends on frequency according to Rayleigh's 
fe 4 -law, so we must replace the Planck spectrum by u(x) = ^L- . 



Figure 1. Decoherence functions (1 — P) for different spectra v and angular 
distributions f(9, ip). The functions have been scaled so as to all have second 
derivative unity at the origin. 



Let us define the coherence time r := A -1 . GRW have given an exact 
solution of Eq. (4) in [1] for the case of a free particle, so to get a feeling 
for what happens we will only mention the simple case when r is is much 
shorter than the dynamical timescale, so that we can neglect the ordinary 
Schrodinger evolution for short times. Thus setting H « and using the 
Gaussian approximation that is valid for t > r, Eq. (4) has the short-time 
solution 



Thus we see that far from the diagonal, for |y — x| » A e //, the elements of 
the density matrix are damped out as e _A * independently of \ e ff- Near the 
diagonal, on the other hand, for |y — x| <C A e //, the crucial parameter is the 




diffusion parameter A = A/A^ 



e//' 



since the damping goes as e 



-A|y-x| 2 t/2 



Table 1 gives A e / / , cj> and r for a variety of radiation sources and Table 
2 gives the diffusion parameter A for the center-of-mass of three different 
objects, in order of decreasing strength. 



Table 1. Properties of various scattering processes 



Cause of collapse 


X eff 




<p[cm. s J 


^electron 


300K air at 1 atm pressure 


0.1 A 




10 24 


10" 13 s 


300K air in lab vacuum 


0.1 A 




10 11 


Is 


Sunlight on earth 


900 nm 




10 17 


6 months 


300K photons 


0.02 mm 




10 19 


1 day 


Background radioactivity 


10" 14 m 


io m 


10" 4 


iu yrs 


I IllflTltllTTl (TYQ^Tlfv 
LlOill L Lllll LlI d V 1 b y 


1 km - 10 


1Q 109 


30s 


CRW pffect 


100 nm 




n/a 


1 0^"\^rt; 

±\J JlD 


r^rtmnir inirrrtwavp Viarkcrnnnrl 


2 mm 




10 13 


10 4 yrs 


Slolar TiPiitrmrm 


0.1 A 




10 11 


10 26 yrs 




3 mm 




10 13 


1 44 vrs 

±\J JlD 


Table 2. Decoherence rate A in 


cm 2 s 1 for various objects 




and scattering processes 










Cause of apparent 


Free 


Wfim Bowling 




wave function collapse 


electron 


dust 


ball 




300K air at 1 atm pressure 


10 31 


10 37 


10 45 




300K air in lab vacuum 


10 18 


10 23 


10 31 




Sunlight on earth 


10 1 


1Q 20 


10 28 




300K photons 


10° 


10 19 


10 27 




Background radioactivity 


10" 4 


10 15 


10 23 




Quantum gravity 


1Q -25 


10 10 


10 22 




GRW effect 


lO" 7 


10 9 


10 21 




Cosmic microwave background 


10 -io 


10 6 


10 17 




Solar neutrinos 


lO" 15 


10 1 


10 13 





The decoherence rates for photons and air agree well with those given 
by Joos and Zeh [11]. The effect of air molecules is seen to dominate at 
room temperature not only at atmospheric pressure but also in a laboratory 
vacuum of 10 6 particles/cm 3 . The radioactivity figures are quite crude, 
since the fluxes of a, (3 and 7 rays vary widely with location and surrounding 
[27]. The energy of the free electron matters only in the case of air [26], 
where it has been taken to be 1 keV. The reason that the cosmic background 
neutrinos [28] are completely negligible compared to those from the sun [29] 
despite similar fluxes is that the weak scattering cross section falls off as 
energy squared in our regime of interest, which is well below the W mass 
of 81 GeV. The neutrino effect is completely impossible to shield against - 
a typical neutrino coasts undisturbed straight through our planet. 

Two sets of numbers from other sources have been put in for compar- 
ison. The GRW values were chosen ad hoc in [1] to match observation as 
well as possible. The quantum gravity values from [30,31] are based on 
dimensional analysis and several perhaps questionable assumptions. For 
instance, A e // is assumed to depend on the mass of the particle, ranging 
from about 1 km for a proton up to astronomical 10 10 m for an electron. 



The main observation to make about the two latter effects is that al- 
though they have been highly publicized, they rank only fifth and sixth 
in strength, trailing by more than twenty orders of magnitude. Hence an 
experiment devised to measure them would have to be nearly perfectly 
shielded from all the stronger sources of decoherence. Blocking out optical 
photons should pose no problem. Background radioactivity could conceiv- 
ably be controlled by using ultra-pure equipment and performing exper- 
iments in a deep mine, thereby avoiding virtually all cosmic rays except 
muons. The effect of air (or any other surrounding substance) and black- 
body radiation from the surrounding is strongly temperature dependent 
(typically A oc T 5 ), and can hence be reduced by nine orders of magnitude 
by working at liquid Helium temperatures. Even under such conditions, 
they would still be stronger than the GRW and quantum gravity effects. 
Although it may become feasible to cool a macroscopically large apparatus 
to microkelvin temperatures, environmentally induced decoherence is in a 
sence endemic: In order to observe the object of an experiment, we must 
by definition let it interact with something else. 

Of no small importance is that we know that all effects in the table 
except those of quantum gravity and GRW do in fact occur. Because of 
this, a number of experiments that have been proposed are likely to yield 
inconclusive results. For example, Hawking [31] and others have conjectured 
that quantum gravity effects might be able to explain the apparent collapse 
of the wave function. Ellis-Mohanty-Nanopoulos have made the estimates of 
such wormhole effects quoted above [30] , but apart from the fact that there 
is no experimentally tested theory of quantum gravity, these effects would 
probably be impossible to detect if an attempt where made to measure 
them, since they would so to speak drown in environmental noise. In fact, 
calculations in [11] show that even ordinary Newtonian gravity often has a 
stronger decoherence effect than quantum gravity. 

The same goes for suggestions to measure an independent GRW effect. 
Squires [21] suggests that a GRW effect might indeed exist and be caused 
by some yet unknown physics that he speculates might be "the physics of 
the 21st century". Our results have shown that the physics of this century 
produces an almost identical effect, and that an additional "new physics" 
contribution of the magnitude postulated by GRW would probably be too 
many orders of magnitude weaker than the decoherence effect to be de- 
tectable. 

Also Rae [22] suggests that the GRW effect might be caused by new 
physics, and speculates that A and a might be new constants of nature. 
Our results have shown that when decoherence dominates, the A and A 
that would be measured in a GRW experiment would be calculable from 
scattering cross sections and spectra of scattering particles. 

The fact that decoherence and the GRW effect have an almost indistin- 
guishable impact on the density matrix opens up an interesting possibility: 
if a GRW effect due to new physics indeed exists, then it might be much 
stronger than originally postulated without contradicting our experience. 
Thus experiments devised to search for a GRW effect are by no means 



without interest. 

4.2. Coherence Lengths 

What is the width Ax of the wavepacket for a free particle? Any 
textbook will give the answer that Ax — > oo as t — > oo, but we are now in 
a position to give a more subtle and indeed finite answer. 

In the presence of our scattering, Ap —> oo like hy/Dt just as for 
Brownian motion in momentum space. For a truly free particle, Ap remains 
constant whereas Ax — ► oo like tAv = tAp/m. An often overlooked fact 
is that there is nothing "quantum" about this whatsoever, other than that 
the uncertainty principle prohibits Ap = initially. Since the free particle 
Hamiltonian is quadratic, the Wigner equation (6) reduces to the Liouville 
equation and the increase in Ax only reflects a classical type of uncertainty, 
our ignorance. 

Let us keep the conventional definition 



for pure states, but redefine Ax for mixed states to be the coherence length, 
roughly speaking the largest distance from the diagonal where the spatial 
density matrix has non-negligible components. Let us redefine Ap analo- 
gously. 

More formally, let us consider a density matrix p that is an incoherent 
mixture of tiltedly squeezed states (Gaussians in Wigner phase space), all 
of which have the same values of the Ax and Ap but with different (x), (p) 
and (xp). GRW analyze the solution to Eq. (1) in detail in [1], and show 
that such a p can be expanded as such a mixture for all times, but with Ax, 
Ap and the mixing function changing with time. During the undisturbed 
Schrodinger time evolution between scatterings, Ax increases as usual while 
Ap remains constant (a shearing in phase space). A scattering causes Ax to 
decrease abruptly while Ap can either increase or decrease. The net result 
of the interplay between these two counteracting effects is that although 
the conventional uncertainties become infinite, both Ax and Ap converge 
to finite limiting values as t — > oo. We might interpret this as that after a 
long time our particle definitely is in a state with spreads Ax and Ap, but 
we have absolutely no knowledge as to where in phase space this state is 
centered. 

In our notation, GRW show that the limiting values of the spreads are 



Ax := (x 2 ) - (x) 



Ax = u \ e ff 



and Ap = v 




T 



where 




and u is the positive solution to 



4u 3 = 2r] 2 uVl + 2u + ?? 4 (1 + u)(l + 2u). 

The dimensionless constant r\ ranges from 10 -29 to 10 26 for the examples in 
the tables above. GRW only consider the macroscopic limit rj <C 1, where 

Ax«2- 1 / 4 ^y /4 and Ap ^ 2^ (h 3 mA) 1/4 , 
but we will also be interested in the microscopic limit r/ S> 1, where we get 
A^2" 1 / 2 4 ! - and Ap^2~ 1 / 2 - h 



mX ef f X ef f 

Table 3 contains the limiting value of Ax for the previously discussed objects 
and decoherence sources. (These figures are to be interpreted as upper 
limits to the true coherence length, since in reality the different effects all 
contribute separately.) "n/a" has been entered for the cases where r is 
greater than the age of the universe, so that the system has not yet had 
time to approach the limiting value of Ax. 

It is interesting to compare the uncertainty products AxAp from above 
with the minimum value h/2 allowed by the uncertainty principle. For the 
macroscopic case i] <C 1, we get 

AxAp r- 

-h/2-~^ 
whereas the microscopic case rj S> 1 yields 



AxAp 
h/2 



~ r? 
^ if . 



This indicates that we should expect to find macroscopic systems such as 
dust particles in states that are nearly minimum uncertainty states, but 
microscopic systems in highly tiltedly squeezed states where the uncertainty 
product is much larger than its minimum. For a free electron decohered 
only by 300K photons, for instance, AxAp « lO 11 ^. This conclusion is 
likely to be valid quite generally, also when the initial state is not of the 
form assumed above, since recent work by Zurek, Habib and Paz [32] has 
indicated that quite general states tend to approach generalized coherent 
states (states with a Gaussian Wigner function) when they interact with 
their environment. 



Table 3. Coherence lengths Ax caused by various decoherence sources 



Cause of apparent Free 10 fim Bowling 

wave function collapse electron dust ball 



300K air at 1 atm pressure 


10" 6 m 


10" 


■17 m 


10" 


-21 


m 


300K air in lab vacuum 


10 7 m 


10" 


13 m 


10" 


-18 


m 


Sunlight on earth 


10 9 m 


io- 


12 m 


10" 


-17 


m 


300K photons 


10 4 m 


io- 


■12 m 


10" 


-16 


m 


Background radioactivity 


n/a 


io- 


11 m 


10" 


-15 


m 


Quantum gravity 


10 4 m 


io- 


9 m 


10" 


-15 


m 


GRW effect 


10 19 m 


io- 


y m 


10" 


-15 


m 


Cosmic microwave background 


10 10 m 


io- 


8 m 


10" 


-14 


m 


Solar neutrinos 


n/a 


n/a 


10" 


-13 


m 



This paper has focused entirely on scattering. The case of coupled 
harmonic oscillators has been studied by numerous authors, and [11] has 
analyzed interactions of Coulomb type. However, many other interesting 
sources of decoherence still remain to be analyzed in detail, and might well 
turn out to be stronger than any of those discussed above. 



5. DECOHERENCE AND THE INTERPRETATION OF 
QUANTUM MECHANICS 

In a discussion of density matrices [33], Feynman writes: "When we 
solve a quantum-mechanical problem, what we really do is divide the uni- 
verse into two parts - the system in which we are interested and the rest of 
the universe. We then usually act as if the system in which we are interested 
comprised the entire universe." In this spirit we summarize our scattering 
results: The effects of "the rest of the universe" can be incorporated into 
two additional terms in the Heisenberg equation of motion; a radiation- 
pressure term that can be absorbed into the Hamiltonian and a dissipative 
term that transforms pure states into mixed states. For macroscopic sys- 
tems, the former is usually negligible whereas the latter effectively damps 
out spatial superpositions. 

As mentioned in the introduction, the decoherence effect removes a 
serious deficiency [34] from non-collapse theories. This makes possible a 
self-consistent interpretation of quantum mechanics that might be called 
the many decohering worlds interpretation: 

Q: Why do we experience collapse? 

A: Because as shown in [3], collapse and QM probabilities will be ex- 
perienced by almost all observers in the grand superposition. 

Q: Why do we only experience macroscopically "nice" superpositions? 

A: Because all others get damped out by decoherence effects before we 
have time to observe them. 

Q: Does the Von Neumann reduction have anything to do with mental 



processes and our observing the system? 
A: No, it has already occurred by the time we observe the object. 
Q: What about Schrodinger's cat? 

A: Due to decoherence, it is for all practical purposes either dead or alive 
- we just don't know which. (Several of the quoted decoherence 
authors would disagree with these answers.) 

If there is agreement on any philosophical interpretation of our cen- 
tury's developments in physics, it is that the ultimate reality is more bizarre 
than anything we ever dreamed of. Hence we no longer feel that we can 
reject a theory merely on the grounds that it involves bizarre and counter- 
intuitive notions. Nonetheless, most people would probably agree that in 
order for any theory to be hailed as fundamental, it should satisfy the fol- 
lowing minimum requirement: it should be able to explain why we do not 
experience any of the bizarre notions that it introduces. 

In the case of quantum mechanics, this would entail using first prin- 
ciples to explain why the wave function seems to collapse and why the 
macroscopical world seems to obey the laws of classical physics. Despite 
Bohr's many examples of classical correspondence for high-energy eigen- 
states, despite Schrodinger's invention of coherent states and despite the 
well-known classical limit when h — > 0, this minimum requirement was not 
met until the work on decoherence in the last two decades. Until then, no- 
body had explained why we never experience say spatial superpositions of 
macroscopic objects (without invoking new physics or superselection rules 
postulated ad hoc). 

While recognizing this success, it is important to remember that al- 
though quantum mechanics including the decoherence effect meets our min- 
imum requirement, this is also all it does. Although the existence of en- 
vironmentally induced decoherence does explain why we never experience 
bizarre macrosuperpositions, and thus makes dynamic reduction mecha- 
nisms (DRMs) such as those proposed by GRW and [18, 35, 36] unnecessary 
for that purpose, it does no more than that. In other words, despite de- 
coherence, macrosuperpositions still do exist - decoherence merely explains 
why we cannot experience them. Instead of being destroyed, superpositions 
spread to ever larger subsystems of the universe as everything gets more 
and more quantum-mechanically entangled with everything else. 

The unease with which many authors view this state of affairs was one 
of the main motivations in introducing DRMs. In DRMs, this perpetual 
entanglement is avoided by making the additional postulate that the wave 
function itself really does get localized, usually in discrete, collapse-like 
jumps. The price is that DRMs need to postulate new physics without 
explaining any phenomena that decoherence alone does not, thus leaving 
themselves vulnerable to Occam's razor. 

As Dieks [37] and others have pointed out, there is also a second source 
of unease. This is that a density matrix describing a mixed state can be 
expressed as a statistical mixture of pure states in infinitely many different 
ways. Thus we are not justified to make the interpretation that an object 



"really" is at a definite position, even if decoherence has made its spatial 
density matrix diagonal. 

But then again, tell one of your friends that the world is a weird place, 
and the answer will be: "So what else is new?" 
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